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Abstract 

When studying the motion of optically trapped particles on the /xs time scale, in low viscous 
media such as air, inertia cannot be neglected. Resolution of unusual and interesting behaviour 
not seen in colloidal trapping experiments is possible. In attempt to explain the phenomena we 
use power spectral methods to perform a parameter study of the Brownian motion of optically 
trapped liquid aerosol droplets concentrated around the critically damped regime. We present 
evidence that the system is suitably described by a simple harmonic oscillator model which must 
include a description of Faxen's correction, but not necessarily frequency dependent hydrodynamic 
corrections to Stokes' law. We also provide results describing how the system behaves under 
several variables and discuss the difficulty in decoupling the parameters responsible for the observed 
behaviour. We show that due to the relatively low dynamic viscosity and high trap stiffness it is 
easy to transfer between over- and under-damped motion by experimentally altering either trap 
stiffness or damping. Our results suggest stable aerosol trapping may be achieved in under-damped 
conditions, but the onset of deleterious optical forces at high trapping powers prevents the probing 
of the upper stability limits due to Brownian motion. 
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I. INTRODUCTION 



In a number of recent publications our group [l 3 1 has presented an experimental system, 
based on the optical trapping of aerosols, for studying Brownian dynamics in both over- 
and under-damped conditions. With this unique perspective we are able to investigate 
the emergence of phenomena such as oscillatory motions due to the influence of inertial 
forces. This is a significant departure from traditional optical trapping experiments which 
are performed under conditions of heavy viscous damping [4|. 

Its understanding |5| can be used in conjunction with optical detection to determine the 

n Q n 

size of microscopic colloidal suspensions [6|,l7| or to measure Avogadro's number [8|. Recently 
it has allowed optical traps to provide a powerful tool in diverse research fields capable of 
acting as a force transducer for molecular biology [9|, viscometry 10J, microscopy [111, and 



fundamental ph ysic s [121 ] . These applications often use the power spectrum method [13j to 



detect position 



14| , measure force s |l5 |. or investigate colloidal dynamics [l6| and rely on 



the study of over-damped systems [17]. 

The investigation of over-damped systems via optical trapping has produced classic exper- 
iments with important ph ysic al results including tests of Kramer's theory 



of critical Casimir forces 



1811, m easurement s 



121 ] and demonstrations of fluctuation theorems 19J. The various 



optical potentials created through optical manipulation have, for example, been used to 



investigate colloidal crystals 



ratchets [22j and freezing 23 1. 



20 



2l|, with particle dynamics providing analogies in thermal 



All experiments in optical tweezers in a liquid environment behave as over-damped os- 



cillators but there have been discussions that under-damped motions are observed 
comments that this should not be possible 



over-damped systems are rare [2 



cle experiments of Ashkin 




17 



2j and 



25| . What is true is that studies of non 



26]. A recent resurgence in the original airborne parti- 
29j looks set to change this, however, with opportunities 
to now study the little probed under-damped systems. The main applications of such ex- 
periments are in aerosol science with most studies to date investigating the chemistry of 



liquid droplets in gas phase environments [30|, |3l| but they also offer more esoteric possi- 



bilities with opportunities to study quantum mechanical effects linked to Brownian motion 



321 ] . Due to the importance of inertia in such systems they provide a drastically different 



damping environment and hence experimental possibilities. 
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When comparing aerosol trapping to aqueous trapping several phenomena may be ob- 
served that would be considered unusual. Varying trapping power alters the axial equilibrium 
position of droplets resulting in 'power gradients' 33] with further increases causing their 
loss. Given a polydisperse nebulised sample the initial power used to capture a droplet 
has pronounced size selectivity [34|, [35[ and once trapped the droplet can undergo vertical 
oscillations at frequencies of ~ 0.1 — 10 Hz. 

This parameter study of the 'mechanical' forces in airborne optical traps was conducted 
in an attempt to explain various phenomena but also proves useful in developing a deeper 
understanding of aerosol trapping. 

A limited exploration has previously been carried out showing a trapped droplet can 



behave in either an over- or under-damped manner and that parametric resonance is easi 
excited In addition we have observed underdamped behavior in two particle systems 



26). 



y 



This paper investigates the transition from over- to under-damped in more detail, exploring 
the parameter space by discussing dependence on laser power, droplet size, and depth into 
the sample at which the droplet is trapped. We show the system can be described by a 
simple harmonic oscillator model by including appropriate corrections. We also test the 
hypothesis that droplets are lost from their traps as the power is increased because they 
cross into the under-damped regime. 

Several methods have been presented to characterise the Brownian dynamics of a trapped 
object, all relying on position measurement, and include, the drag force method, the 
equipartition method, the step response method, autocorrelation, and the power spectrum 
method [36j. Position sensitive detectors, video tracking or quandrant photodiodes (QPDs) 
are often used to detect particle position each with their own speed and precision advantages 



and disadvantages 
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381 ] . Here we employ a QPD due to their high bandwidth and use 



the power spectrum method, considered the most reliable |39|], to characterise the Brownian 
motion of droplets within optical traps. The autocorrelation of a single particle is also feasi- 
ble but can provide poor results in noisy systems. The power spectrum method decomposes 
the motion into frequency components so any noise can be easily dealt with. Although we 
are mainly concerned with observing the dynamics of trapped aerosols it is feasible to use 
this method to measure precise forces and position. Normally the method allows calculation 
of trap stiffness with prior knowledge of the viscosity of the surrounding medium and the 



particle radius 36], but it will be shown that when in air only the radius is needed. 
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First we will discuss the theory used to describe our experimental system along with any 
subtle corrections that may need to be considered and what the magnitude of their effect 
would be. We will then describe our experimental apparatus and procedures which differ 
slightly from optical trapping in aqueous media. We will present evidence that supports 
that our system is described by a simple harmonic oscillator model with certain corrections 
necessarily included. Finally we will try to determine the cause of the unusual behaviour 
partly by studying how trapped droplets behave near the critically damped regime. 

II. THEORY 

Throughout this work we assume the particle velocity is well below the speed of sound and 
the propagations of interactions in the fluid are instantaneous, hence the fluid, air, is treated 
as incompressible 40(. An optically trapped particle is treated as residing in a harmonic 
potential well experiencing a Hookean restoring force when displaced through Brownian 
stochastic forces. The characteristic time for such a particle to lose energy through friction 



is a balance between inertial and viscous forces, ti nert = m/70 |41|, where m is the particle 
mass and 70 is the viscous drag. For the smallest droplet studied here this time is longer 
than our experimental resolution so inertia must not be neglected as is usual for studies in 
liquid media. The Langevin equation describing the motion of a liquid aerosol of radius R, 
mass m, optically trapped in a fluid of temperature T, kinematic viscosity u, and density 



P fluid, with stiffness k is [42] 



x(t) + Tx(t) +tt 2 x(t) = A({t), (1) 



where f2 = y(ft/m) is the natural angular frequency of the droplet position fluctua- 
tions, T = 6irr]R/mCr is the viscous damping of the medium due to a viscoscity 77, 
A = (2fc#Tr /m) 1 / 2 43j, |44| is the Brownian stochastic force where ks is the Boltzmann 



constant, and for all t and t' (((t)) = and (((t)((t')) = 5(t — t'). Stokes' Law is corrected 
for finite Knudsen number effects by including the empirical slip correction factor, C c , with 



a 5.5% — 1.6% reduction in drag for 3 — 10/xm diameter droplets, respectively [45j. Fourier 
transforming equation [1] and finding the expectation value we decompose the motion into 
frequency components and find the power spectrum of position fluctuations to be 
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where u; is the angular frequency. This spectrum has a characteristic high frequency tail 
with bj~ A gradient and a plateau value at low frequencies equal to 2ksTY / kQ 2 . As inertia 
is included there is an additional limiting case, compared to over-damped oscillators, at the 
point of inflection equal to 2kBT/nT. We define the ratio of damping coefficient to natural 
frequency as the 'damping ratio', T/Q. For over-damped systems this is always greater than 
unity, as is found for colloidal systems where it is usually greater than ten. In such cases 
the first angular frequency term in the denominator can be neglected with respect to O to 
give the usual power spectrum for over-damped optical traps [2] . When trapping in air the 
system has the potential to become under-damped and hence T/Q < 1. 

The above Langevin equation assumes the motion occurs in bulk fluid media with uniform 
velocity, far away from other objects and surfaces so Stokes' law is only corrected for finite 
Knudsen number. However, in reality the objects here are undergoing linear harmonic 
motion within significant proximity (< 10 R) of a coverslip so it is inappropriate to assume 
Stokes' law still applies. Studying this problem it is seen there are two significant corrections 
that may need to be applied to the 'in bulk' theory. 



Firstly, the object is undergoing linear harmonic motion and so a more comp 



dependent, friction must be considered which was shown by Stokes to be [40 



ex, frequency 



41 



46| 
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The first term comprises the familiar Stokes' drag plus a frequency dependent correction. 
The second term arises from the inertia created by any fluid entrained due to the past motion 



481 ] . at times 



of the particle. This hydrodynamic correction is often neglected j^, [9, [m| 4^, 
with good cause, but needs to be applied when requiring precision > 10% [4l|. Here we will 
try to justify our exclusion of such terms in a little detail. 

Using equation [3] and following Berg-S0rensen and Flyvbjerg J4I] we derive the hydrody- 
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namically correct power spectrum in angular frequency to be 
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where u m = T/(l + 2 ' Kp ^ uidR ) an d u v = 2v/R 2 . For systems incorporating inertia the usual 
definition of a corner frequency, o; c = re/70, clearly can no longer apply so the dependence 
has been removed. The low density of air reduces the denominator of u m to close to unity, 



effectively removing any effective mass considerations as the entrained fluid is negligible 



49]. 



Decomposing equation [3] into frequency components via Fourier theory (equation 31 in Berg- 



S0rensen and Flyvbjerg 



411 ]) shows the larger kinematic viscosity of air, hence u u , reduces 



any correction to Stokes' law compared to trapping in water. 

The effect on the hydrodynamic correction of a negligible effective mass and larger u u 
is not immediately apparent. To visualise the relative size of the correction we plot fig. [U 
the ratio of the power spectrum in equation [2] to the hydrodynamically correct version in 
equation HI 




10 100 1000 

Angular Frequency/2ji (Hz) 



10000 



FIG. 1: Theoretical plot of S X J S x - hydro as a function of angular frequency. For a given particle 
radius the hydrodynamic correction is smaller in air (red dot dashed) than in water (blue dashed), 
but, for the large liquid aerosols (black solid) the effect starts to become significant. Trap stiffness, 
k = 2 pN^m -1 . 

Clearly for a given particle type and size the correction is significantly smaller when 
studying aerosols. However, the solid line shows the error begins to become significant 
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for aerosols with a radius that would be considered relatively large for particles normally 
used in power spectrum based studies in liquid. Should further studies be performed the 
hydrodynamic correction must be investigated to improve accuracy and precision. 

For the majority of applications optical trapping parameters, such as trap stiffness, are 
only needed with an accuracy of ~ 10%, and considering the magnitude of the correction 
factor we believe it is reasonable to neglect the effect of frequency dependent friction in this 
study. 

The second correction to be considered is that given by Faxen regarding the force on a 
sphere in motion near a plane surface, exactly what occurs when trapping with high NA 
optical tweezers due to the proximity of coverslips. Here we only consider the correction in 
the lateral direction although both axial and rotational equivalents exist 
law shows the viscous drag on a sphere increases as it approaches a plane surface according 
to [52J 



50 



51] . Faxen 's 



1 _ ( 1K\ _J_ I ( R\ _ JL ( Ry 1_ ( R\ 

1 {wLJ^syLJ 256 \L) 16 \LJ 

where L is the distance between sphere centre and surface. For the particle sizes studied here 
this can have a dramatic effect on the friction experienced; even when trapping at distances 
approaching 40 /im from coverslips there can be a 7% increase. 

The final theoretical consideration is that to compare power spectra in given data sets 
we must calculate the detection system sensitivity, ft, given in volts output per unit dis- 
placement of the particle. This is because the sensitivty can alter between experiments due 
to variations in power or simply geometry at the focus. Finding ft allows voltage power 
spectra, those recorded directly from the experiment, to be converted to physical spectra, 
nm 2 Hz _1 versus Hz. Conventional methods rely on the relative simplicity of colloidal sys- 
tems by using, for example, the drag force method [53] , its extension to an oscillating samp 



stage 



54j or moving a fixed bead over a known distance through the laser beam waist 55]. 



Clearly the former two would be difficult to implement in air and the latter is obviously not 
a good replica of experimental conditions [ssj] . A recent technique has been demonstrated 
that combines two methods to measure detector calibration from experimentally measured 
values alone j^] . It is hoped, even with the unique problems of airborne trapping, by using 
AODs or SLMs to oscillate the trap position this technique will be developed for future 
experiments. 
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Here we are not concerned with high precision and for simplicity we calculate the detector 
sensitivity, (3, from an uncalibrated voltage power spectrum S v (oj) = f3 2 Si nert {<jj) using the 
plateau value, P v , reached for us ^> Q in the function u 4 S v (u) |l3[. We find the detector 
sensitivity, 0, to be 



P v m 



III. EXPERIMENTAL 



Droplets are trapped using a custom built inverted tweezers pictured in fig. |2j The 
beam from a 532nm Laser Quantum Finesse 4W c.w. laser is expanded by a Keplerian 



telescope to slightly overfill 57] the back aperture of a Nikon Plan lOOx (NA = 1.25 58]) 
oil immersion microscope objective. The beam is focused though a type one cover slip 
into an aerosol chamber constructed from a cylindrical plastic enclosure 9mm in height and 
35mm in diameter. This produces an enclosed environment where a high relative humidity 
can exist and also shields the trapping region from external air currents. The top of the 
chamber is made from a type zero cover slip to allow for transmission and then collection of 
the scattered trapping laser by a long working distance (LWD) Mitutoyo lOOx (NA = 0.55) 
objective, whose back aperture is imaged equally onto the four quadrants of a quadrant 
photodiode (QPD) (Hamamatsu Silicon Diode Array S5980) via a 4f imaging system. The 
Mitutoyo objective also acts as the condenser lens for Kohler illumination (not shown). The 
Nikon objective and an appropriate tube lens images the sample through a laser filter onto 
a Basler A602f firewire camera. 

The liquid aerosol is produced by nebulising a salt solution (20-80g/L) with an Omron 
MicroAir NE-U22 vibrating mesh nebuliser which produces a polydisperse sample of liquid 
droplets with a mass median aerodynamic diameter of 4.9/xm [59]. The aerosol is transferred 
through a hole in the chamber side via a custom made tapered glass nozzle [so | . 

The trapping beam is focused ~ 30/im above the coverslips which are soaked in a 50% 
aqueous dilution of 'Decon 90' for longer than one week. This treatment increases the 
hydrophilicity of the glass and once aerosol has deposited on the cover slip it provides a 
relatively thin, flat, and uniform film of water above which we trap. De-ionised water 
saturated tissue paper is also placed in the chamber to increase the relative humidity, but 
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FIG. 2: Apparatus diagram. A Gaussian beam is expanded by lenses LI and L2 and directed to 
slightly overfill the back aperture of the Nikon objective (TOBJ) with mirrors M and DM. The 
long working distance objective (LWD) collects the scattered light from the droplet and its back 
aperture is imaged onto the QPD via a 4f lens system. Power is controlled using a polarising 
beam cube (PBC) and half wave plate (WP). The same Nikon objective with an appropriate tube 
lens (TL) is used to image the sample (S) through a dichroic mirror (DM) and filter (F) onto the 
firewire camera (CMOS). The QPD, COB J, and TOBJ are each mounted on three axis translation 
stages with the axial axis of TOBJ controlled either manually or by digital micrometer. BD is a 
beam dump. 



we ensure it does not touch the cover slip as this can induce flows in the water layer. 
Figure |3] shows an enlarged view of the trapping region geometry and also explains the 
relation between trapping height, L, and objective displacement, X. 

Control over droplet size was required to fully investigate observed phenomena. Firstly, 
this was achieved imprecisely by varying the concentration of the nebulised salt solution [61 1 
as a higher concentration decreases the droplets vapour pressure allowing them to equilibrate 
with their surroundings at larger sizes. Secondly, more precise size selectivity can be induced 
with, on average, a positive linear dependence of captured droplet size on laser power 3^, 35 ]. 

Having trapped a droplet the nebuliser is turned off. Once the droplet has reached 
equilibrium with its surrounding environment, and the remaining aerosol settled, the current 
produced by the detection of light on the QPD is sent, via shielded cables, to amplification 



electronics 



62] containing a 50 kHz anti-aliasing filter. Data was acquired at a sampling 



frequency of 50 kHz for four seconds with a National Instruments PCI-6014E DAQ card, 



9 











L=(X^-Ah) r ^ L 
n g n w 


L 


Air 


n 3 


, lAh 


Aqueous 


n w 




Glass 


n g 





FIG. 3: Enlarged view of the trapping region in the sample of figure [2l The refractive indices of the 
coverslip and index matched oil, water, and air are n g , n w and n a respectively. Ah is the thickness 
of the water layer. Displacing the objective X microns from being focussed on the first interface 
displaces the particle a distance L, given in the figure. 



in differential mode. The voltage difference between left and right pairs of quadrants on 
the QPD represents the x position and the difference between the top and bottom pairs 
represents the y position. The voltage versus time data was Fourier transformed using 
Lab VIEW and all remaining data analysis was performed offline at a later time. In order 
to minimise any parameter variation over time the experiments were carried out as quickly 
as possible with raw voltage versus time data not saved to increase speed still further. The 
detailed analysis of the data obtained for a colloidal case is extensively described in Berg- 
S0rensen and Flyvbjerg ^1|, and much remains the same here. An image of the trapped 
droplet was also taken with each power spectrum for later analysis. 

To further reduce background noise work was always carried out, where possible, solitarily 
in the laboratory with the laser used at > 30% capacity and power control achieved by 
using a pair of half wave plates with polarizing beam cubes. The first split the beam for 
two different experiments and the second controlled power for this experiment alone. The 
power was varied between a minimum of 0.702 ± 0.009mW and a maximum of 510 ± 6mW. 

Unlike tweezing in water, simply increasing the trapping power does not assist in cap- 
turing an aerosol droplet from the nebulised cloud and as such initial laser power must be 
carefully selected. Each droplet trapped was subject to an increase in laser power in uniform 
steps with power spectra measurements taken at each. The minimum attainable damping 
ratio for each droplet was taken from the last power spectra measured before it fell from the 
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trap upon increasing the power (i.e. the highest power). This represents an upper limit on 
the ratio for that size. 

To study how the water-air interface to droplet height may affect the dynamics we simply 
keep a constant laser power and vary the height of the sample stage, controlled and measured 
by micrometer. The water layer thickness was measured by observing when a reflection 
of the trapping beam focus is obtained at both the water-air and glass-water interfaces. 
Having been focussed through two refractive index mismatched interfaces (glass to water 
and water to air) there will be an associated focal shift 63[ of which a rigorous description 



is complex 



64. 



65| and not discussed here. A simple paraxial approximation is used to 



calculate the droplets position inside the chamber given a vertical displacement of the sample 
stage around a fixed objective. Modelling of the axial equilibrium position of the trap 
and experiments imaging the droplet from the side indicate that the relationship between 
droplet height and objective displacement is linear (data not shown) supporting the paraxial 
assumption |33l. 166]. 



IV. RESULTS AND ANALYSIS 



Typical power spectra of position fluctuations from optically trapped droplets are shown 
in figure HI illustrating, for a 3.7 ± 0.2 /xm radius droplet, the ease with which the system 
can be transferred between over- and under-damped dynamics by varying laser power. The 
tail falls off with oj~ 4 as expected for u Q from equation [2] and a clear resonance peak 
begins to establish itself, indicative of the droplet moving through the critical and into the 
under-damped regime. 

For completeness, a plot of the autocorrelation function [16( of a single droplet in an 
under- and over-damped state is shown in figure [5j It shows the classic exponential decay 
for over-damped motion and sinusoidal oscillation enveloped by exponential decay for under- 
damped oscillators as expected. 

The trend seen in figure H] remains for all droplets; an increase in power increases lateral 
trap stiffness and moves the system towards or into the under-damped regime. A range of 
damping ratios has been observed from 3.57±0.07 down to 0.260±0.006 over the 4.7±0.5 /im 
radius range studied. There is also an associated decrease in area under the power spectrum 
curve with increasing laser power, indicating a reduction in the position variance of the 
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FIG. 4: Power spectra of a droplet of radius 3. 7±0. 2 trapped at powers 40. 9±0. 5 (circles), 130±2 
(squares), and 356 ± 4 (triangles) mW resulting in damping ratios of 1.69 ± 0.04, 0.794 ± 0.01, and 
0.364 ± 0.001 respectively. The fitting parameters for the top (blue), middle (purple) and bottom 
(red) curves are 0/2vr = 690±13, 1565±7, 3106±13 Hz and r/2vr = 1167±13, 1242±17, 1131±13 Hz 
respectively. As the power increases the appearance of a resonance peak is clear indicating the 
move into an under-damped regime, along with a decrease in area and hence position variance. 
The natural frequency increases with laser power because of the associated increase in lateral trap 
stiffness, n. 

droplet. 

The inclusion of inertial terms in the Brownian theory means only the mass of the particle 
is needed to calculate trap stiffness. Using the radius from video microscopy we obtain lateral 
trap stiffness values ranging from 0.12±0.10 to 98±17 pN/zm -1 for 1.0±0.3 to 5.7±0.4 jum 
radius droplets. 

One would expect the natural frequency of trapped droplets to vary as the square root 
of laser power, assuming the trap stiffness is linearly proportional to trapping power. This 
is confirmed in figure [6] for a 1.8 ± 0.2 /zm radius droplet. For the range of radii and powers 
studied here we observe natural frequencies between 27r(328 ± 12) rads -1 and 27r(3433 ± 
15) rads -1 , falling close to and well above the corner frequencies measured by tweezers in 
liquid based systems, although obviously not directly comparable. 

The above results produce a downward shift in damping ratio by increasing lateral trap 
stiffness with larger laser powers. A decrease in friction felt by the droplet could likewise 
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FIG. 5: Plot of experimental autocorrelation data with associated fits for a 5.2 ± 0.2 /jm optically 
trapped aerosol in an over- (triangles) and under-damped (circles) state trapped with powers 74±1 
and 442 ± 7 mW respectively. Both traces clearly follow the classic exponential decay except in 
the under-damped case there is also the sinusoidal oscillation expected. 
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FIG. 6: An example of how natural frequency for a 1.8 db 0.2 //m radius droplet increases with the 
square root of laser power as expected from f2 = yjn/m. The lateral trap stiffness axis is displayed 
for interest and is non linear. The error bars are standard error of the mean for the natural 
frequency rather than the trap stiffness (although they are smaller than the points themselves). 

shift the ratio by varying the damping and as stated earlier Faxen's correction predicts that 
the proximity of a surface to our microscopic object heavily influences this. Utilising this 
surface to droplet height dependence figure [7] demonstrates that lowering the sample stage, 





- 0.78 



13 



hence increasing the distance, reduces the damping and transfers the system from over- to 
under-damped. Note the resonance peak remains approximately at the same frequency for 
each spectrum as only the damping is changing, not the trap stiffness, contrary to figure HI 
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FIG. 7: Power spectra demonstrating changes in damping as a function of surface-droplet distance 
for a 3.8 ± 0.2 /im radius droplet trapped with 46.3 ± 0.6 mW. The droplet was moved to heights 
of 4± 1 /im (circles), 9± 1 /xm (squares), and 14 ± 1 /im (triangles) above the water layer resulting 
in damping ratios of 3.40 ± 0.06, 1.06 ± 0.01, and 0.92 ± 0.01 respectively. The fitting parameters 
for the top, middle and bottom curves are Q/2ir = 981 ± 5,962 ± 5,815 ± 10 Hz and T/2ir = 
904 ± 10, 1019 ± 12, 2773 ± 32 Hz respectively. The middle and top spectra are multiplied by 25 
and 200 respectively to displace the data on the y-axis for clarity. 

Extracting damping values from data similar to fig. [7J we can plot the dependence of 
friction upon droplet-surface height to obtain fig. [8] [1]. Here the micrometer raised the 
sample stage in increments of 1/im, decreasing to 0.5/^m as the surface was approached. 

In fig. [9] we plot the natural frequency as a function of height from the water layer. There 
is a steady fall off with distance indicating the spherical aberration induced is degrading the 



561 ] the data is approximately linear as 



trap stiffness. Unlike the data of Vermeulen at al. 
the trap stiffness is independent of the viscous damping. 

Figures 0H6] and [9] confirm the simple harmonic oscillator model applies to our experimen- 
tal system for all regimes of damping. Our results are not precise enough to examine the 
need for frequency dependent hydrodynamic correction, but figs. [7] and |8] do show Faxen's 
correction has an important effect on the damping experienced. 
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FIG. 8: Variation of damping experienced by a droplet, trapped with 6.2 ± O.lmW, as a function 
of surface-droplet height. 
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FIG. 9: Natural Frequency as a function of distance from water layer for the droplet in fig. 7. 

Bearing this in mind we turn to the most curious unusual phenomena observed; the loss of 
droplets from traps at a particular upper limit on laser power. The confirmation of a simple 
harmonic oscillator model as an appropriate description of our system leads us to test the 
hypothesis that the particles become largely under-damped so quickly that instability is 
caused. To test this and ascertain the true cause we plot, in figure [I0J the upper limit on 
the minimum damping ratio attainable against droplet radius. 

Figure [HD illustrates droplets can exist in an under-damped regime. Some droplets are lost 
from traps while over-damped and some while under-damped. This suggests no instability is 
induced as the object crosses the critical damping point and so is not the reason for droplet 
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FIG. 10: Upper limit on the minimum attainable damping ratio against droplet radius. It is 
an upper limit as we increase the laser power in finite increments. The dashed horizontal line 
represents a critically damped system. The error bars are standard error of the mean. 

loss with increasing trapping power. There is clear size dependence but one must be careful 
to note this does not lead to the conclusion that droplets do indeed become unstable as 
they reach a particular damping threshold, but rather fall from the optical traps at a certain 
upper limit on trapping power. 



V. DISCUSSION 



The reader may notice significantly more scatter in the results of this investigation com- 
pared to experiments using similar techniques in liquid. The reasons will now be discussed 
briefly showing the complex nature of the experiment and the engineering challenges faced 
to improve future precision. 

The majority of previous work using the power spectrum method is based on tweezing 
solid microspheres, with precisely known radii, in the liquid phase allowing very high preci- 
sion studies; indeed, the ability to detect sphere non- uniformity is possible [54j . In the studies 
here a large source of error is measuring the radius using video microscopy with the likely 
errors propagating heavily into some of the systems calculated properties (k oc m oc R 3 ). 
Trappin g of solid aerosols, with known radii would remove this problem but this is more 
difficult 6oj|. In addition, highly precise radius measurements are possible via CERS 67] 
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but require sensitive spectrometers. We have shown in a previous publication that fitting 
to the data of fig. [8] with equation [5] allows the droplet radius to be measured with good 
precision l|. 

It is difficult to determine which individual factor, trap stiffness or damping, contributes 
to the variation in damping ratio for any given experiment. For an individual droplet the 
surrounding conditions can remain relatively constant over the time of a single experiment 
as, with no additional aerosol flow from the nebuliser, the droplet quickly reaches equilibrium 
with its surroundings. To trap another droplet nebulisation must resume where upon the 
chamber conditions can alter. Additional aerosol can settle on the coverslip changing the 
thickness of the aqueous layer and hence the optical potentia 
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at the trap site 



6611 . As mentioned variation 



331 ] and hence distance 



together with the proximity of the particle to the surface 
in trapping power between droplets alters the height of the droplet 
from the underlying water layer, thus again altering the optical potential and damping. 
These factors contribute to the rather complex and difficult analysis of the system. 

The difficult nature, relative to colloidal tweezers, of trapping in air imposes several im- 
portant experimental methods. A long working distance condenser must be used due to the 
aerosol chamber height, but a higher NA lens may have been desirable to improve detector 
sensitivity [70|. Most colloidal experiments use monodisperse suspensions of solid particles 



thus allowing an arbitrary number of measurement repetitions; often up to 100 power spec- 
tra are averaged. However, with the dynamic system investigated here (the droplets are 
continuously finding an equilibrium with the surrounding environment) the conditions of 
the experiment may not remain constant long enough for repeated measurements to im- 
prove precision, hence the choice of sampling and no averaging over multiple power spectra. 
Also, we are looking at an inherently unstable region with the aim, at times, of losing the 
trapped droplet so, clearly, another particle of the exact same size and composition cannot 
be found. With the current iteration of apparatus there is a clear trade off between speed 
and precision. 

Some studies have used a secondary, independent probe beam to monitor position fluc- 



tuations as this allows greater flexibility and perhaps improved accuracy [71j. We employ 
only a single beam because a very small amount of power is needed to tweeze in air 35) and 
a second beam has the potential to significantly alter the potential at the trap site. 

In future studies we suggest that a system including a 'science chamber' be developed 
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where many variables can be controlled. A particle could be trapped and transferred to 
such a chamber with relative humidity control, with or without a water layer, and with a 
lower physical profile to enable the use of higher NA condenser optics. Also the mechanical 
stability of our system is not fully optimised so the precision seen could be improved. 

Not discussed in detail here is that Faxen derived his correction for a sphere moving with 
constant velocity which is not the case, so for a complete solution the frequency dependent 
friction should be combined with Faxen's correction 4l|, |54( . In addition when dealing with 
solid in fluid systems no slipping occurs at the boundary between the two materials upon 
translation. However, the physics involved becomes more complicated when studying fluid 
in fluid systems; slip can occur. Due to the possibility of slip at the surface of the fluid 
sphere, flow can be induced inside the water droplet. This flow causes reduction of the well 

52! . Q according to 



known pre-factor of Stokes' law 

_ QnvpRv 1 + f <r 

■f stokes 1 ; \l ) 

C c l + a 

where a is the ratio of the dynamic viscosities of the medium and droplet, \x m and fi p 
respectively, giving Stokes' law for a water droplet in air, to be 

5.967rupRv 



stokes 



c c ' 

which has been taken into account in the data analysis here. 



(8) 



VI. CONCLUSION 



The work here is the first parameter exploration of the Brownian motion of optically 
trapped liquid aerosols. We have presented evidence that the system is suitably described by 
a simple harmonic oscillator model which must include a description of Faxen's correction, 
but not necessarily frequency dependent hydrodynamic corrections to Stokes' law. The 
results also show there is difficulty in decoupling the parameters responsible for the observed 
behaviour. Having hyothesised that an instability is caused in the system when crossing from 
over- to under-damped regimes we see this is not supported by the evidence. 

Considering the Langevin equation it is seen that there are only four processes providing 
forces that give rise to droplet position fluctuations; Brownian white noise, friction, inertia, 
and the optical force. Having rejected any damping or inertial cause for the instabilities with 
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the evidence presented here the logical conclusion is that the optical force must determine 
whether the droplet remains trapped or not and gives rise to the size dependence of fig. [TQl 
This will be discussed in a further publication. 

The investigation has provided results extending the boundaries of precise studies of 
Brownian motion in optical tweezers into a new damping regime. It is hoped these results 
will provide researchers with a new understanding of optical tweezers for studies in both 
fundamental and applied science, providing a rich playground of study in the under-damped 
regime. 
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Abstract 

When studying the motion of optically trapped particles on the lis time scale, in low viscous 
media, inertia cannot be neglected. Resolution of unusual and interesting behaviour not seen in 
colloidal trapping experiments is possible. Using power spectral methods we explore, in detail, the 
Brownian motion of optically trapped liquid aerosol droplets concentrated around the critically 
damped regime. We show that due to the relatively low dynamic viscosity and high trap stiffness 
it is easy to transfer between over- and under-damped motion by experimentally altering either 
trap stiffness or damping. The difficulty in decoupling all the parameters responsible for the 
observed behaviour is illustrated and discussed. Furthermore, it is believed these experiments show 
promise for analogies in many areas of physics and contribute significantly to the understanding 
for applications beyond this field. 
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I. INTRODUCTION 



Brownian motion provides many opportunities for analogies in areas beyond microscopic 
diffusion, for example in the stock market [1], and over fOO years after its explanation [2] 
research continues into its further understanding and application [3]. Its understanding can 
be used in conjunction with optical detection to determine the size of microscopic colloidal 
suspensions [4, 5] or to measure Avogadro's number [6]. Recently it has allowed optical traps 
to provide a useful tool in diverse research fields capable of acting as a force transducer for 
molecular biology [7], viscometry [8], microscopy [9], and fundamental physics [10]. These 
applications often use the power spectrum method [11] to detect position [12], measure 
forces [13], or investigate colloidal dynamics [14] and rely on the study of over-damped 
systems [15]. 

The investigation of over-damped systems via optical trapping has produced classic exper- 
iments with important physical results including tests of Kramer's theory [16], measurements 
of critical Casimir forces [10] and cross correlations between colloids [14]. The various optical 
potentials created through optical manipulation have, for example, been used to investigate 
colloidal crystals [17, 18], with particle dynamics providing analogies in thermal ratchets [19] 
and freezing [20]. 

All experiments in optical tweezers in a water type environment behave as over-damped 
oscillators but there have been discussions that under-damped motions are reachable [21] and 
comments that this should not be possible [15, 22]. Rarely have systems that are not over- 
damped been investigated [21, 23] but there is a recent resurgence in the original airborne 
particle experiments of Ashkin [24-26] . Their main application is in aerosol cience with most 
studies to date investigating the chemistry of liquid droplets in gas phase environments [27, 
28]. Due to the importance of inertia in such systems they provide a drastically different 
damping environment and hence experimental possibilities. 

Few investigations into the mechanical dynamics of airborne systems [29] have been per- 
formed but their study is important with the realisation they could provide access to the 
little probed, under-damped regime [23] and new unexploited applications. 

During trapping in air, when compared to trapping in liquid media, several phenomena 
may be observed that would be considered unusual. Varying trapping power alters the axial 
equilibrium position of droplets resulting in 'power gradients' [29] with further increases 
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causing their loss. Given a polydisperse nebulised sample the initial power used to capture 
a droplet has pronounced size selectivity [30, 31] and once trapped the droplet can undergo 
vertical oscillations at frequencies of ~ 0.1 — 10 Hz. 

The aim of this study is to develop an understanding of the process of aerosol trapping 
and to see if the unusual phenomena can be explained. We have previously shown that single 
optically trapped droplets can behave in either an over- or under-damped manner and that 
parametric resonance is easily excited [23]. Here we look at the transition from over- to 
under-damped in more detail, exploring the parameter space by analysing dependence on 
power, size, and depth into the sample. To this end we desire to characterise the dynamics 
of an optically trapped droplet. Several methods have been presented to measure the forces 
exerted on a trapped object, all relying on position measurement, and include, the drag 
force method, the equipartition method, the step response method, autocorrelation, and the 
power spectrum method [32]. 

Position sensitive detectors, video tracking or QPDs are often used to detect particle 
position each with their own speed and precision advantages and disadvantages [33, 34]. 
Here we employ a QPD due to their high bandwidth and use the power spectrum method, 
considered the most reliable [35], to characterise the Brownian motion of droplets within 
optical traps. The autocorrelation of a single particle is also feasible but can provide poor 
results in noisy systems. The power spectrum method decomposes the motion into its 
frequency components and so any noise can be easily dealt with. Although we are mainly 
concerned with observing the dynamics of trapped aerosols it is feasible to use this method to 
measure precise forces and position. Normally the method allows calculation of trap stiffness 
with prior knowledge of the surrounding medium's viscosity and the particle radius [32], but, 
it will be shown that when in air, only the radius is needed. 

The results provide the first tentative steps towards using optical tweezers to probe com- 
plex under-damped systems where inertial effects can no longer be neglected; providing a 
rich playground of study in the under-damped regime. 

II. THEORY 

Throughout this work we assume the particles' velocity is well below the speed of sound 
and the propagations of interactions in the fluid are instantaneous, hence we treat the fluid, 
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air, as incompressible [36] . An optically trapped particle is treated as residing in a harmonic 
potential well experiencing a Hookean restoring force when displaced through Brownian 
stochastic forces. The characteristic time for such a particle to lose energy through friction 
is a balance between inertial and viscous forces, t inert = m/7 [37]. For the smallest droplet 
studied here this time is longer than our experimental resolution so inertia must not be 
neglected as is usual for studies in liquid media. The Langevin equation describing the 
motion of a liquid aerosol of radius R, mass m, optically trapped in a fluid of temperature 
T, kinematic viscosity v, and density p fluid, with stiffness k is [38] 



where Q = y(/t/m) is the natural angular frequency of the droplet's position fluctuations, 
T = 67ir]R/mC c is the viscous damping of the medium, and A = (2/csTr /m) 1 / 2 [39, 40] is 
the Brownian stochastic force where ks is the Boltzmann constant. Stokes' Law is corrected 
for finite Knudsen number effects by including the empirical slip correction factor, C c , with 
a 5.5% — 1.6% reduction in drag for 3 — 10/im diameter droplets, respectively [41]. Fourier 
transforming equation 1 and finding the expectation value we decompose the motion into 
its frequency components and find the power spectrum of position fluctuations to be 



where u is the angular frequency. This spectrum has a characteristic high frequency tail 
with u~ A gradient and a plateau value at low frequencies equal to 2ksTT j 'kVL 2 . As inertia is 
included there is an additional limiting case, over the over-damped oscillator, at the point of 
inflection equal to 2/c#T / kY. We define the ratio of damping coefficient to natural frequency 
as the 'damping ratio', T/Q. For over-damped systems this is always greater than unity, as 
is found for colloidal systems where it is usually greater than ten. When trapping in air the 
system has the potential to become under-damped and hence T/Q < 1. 

The above Langevin equation assumes the motion occurs in bulk fluid media with uniform 
velocity, far away from other objects and surfaces so Stokes' law only needs to be corrected for 
finite Knudsen number. However, in reality the objects here are undergoing linear harmonic 



x(t) + r±(t) + n 2 



x(t)=Ar)(t), 



(1) 
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motion within significant proximity (< 10 R) of a coverslip so it is inappropriate to assume 
Stokes' law still applies. Analysing this problem it is seen there are two significant corrections 
that need to be applied to the 'in bulk' theory. 

Firstly, the object is undergoing linear harmonic motion and so a more complex, frequency 
dependent, friction must be considered which was shown by Stokes to be [36, 37, 42]; 



friction 



-to 1 1 + U - 1*p^R 3 i + 5 Vis I *■ (3) 



The first term comprises the familiar Stokes' drag plus a frequency dependent correction. 
The second term arises from the inertia created by any fluid entrained due to the particle's 
past motion. This hydrodynamic correction is often neglected [7, 21, 23, 43, 44], at times 
with good cause, but needs to be applied when requiring precision > 10% [37]. Here we will 
try to justify our exclusion of such terms in a little detail. 

Using equation 3 and following Berg-S0rensen and Flyvbjerg [37] I derive the hydrody- 
namically correct power spectrum in angular frequency to be 

,1/2N 

•j ! - i ; !i i 

9k, T 

-tinert ■ - ^rvh-L 



n 2 r [i 

ST-hydroi 10 ) = ~ } / N ~\2 / / U 2' (^) 



where oo m = T/(l + 2 *£^ ) and uo v = 2u/R 2 . For systems where inertia is included the 
definition of u c clearly can no longer apply so the dependence on corner frequency has been 
removed. The low density of air reduces u; m 's denominator close to unity, so effectively re- 
moving any effective mass considerations as the entrained fluid is negligible [45] . Decompos- 
ing equation 3 into frequency components via Fourier theory (equation 31 in Berg-S0rensen 
and Flyvbjerg [37]) shows air's larger kinematic viscosity, hence cu u , produces a reduction in 
any correction to Stokes' law compared to trapping in water. Considering these arguments, 
reviewing relevant and similar work, and understanding we are not attempting to make 
highly precise measurements, we believe we are justified in neglecting the frequency depen- 
dent friction correction. Also note that other sources of error inherent in the experiments, 
discussed later, will dominate over any induced by this exclusion. 

The second correction to be considered is that given by Faxen regarding the force on a 
sphere in motion near a plane surface, exactly what occurs when trapping with high NA 
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optical tweezers due to the proximity of coverslips. Here I only consider the correction in 
the lateral direction although both axial and rotational equivalents exist [46, 47]. Faxen's 
law shows the viscous drag on a sphere increases as it approaches a plane surface according 
to [48] 




where L is the distance between sphere centre and surface. For the particle sizes studied here 
this can have a dramatic effect on the friction experienced; even when trapping at distances 
approaching 40 /im from coverslips there can be a 7% increase. 

The final theoretical consideration is that to compare power spectra in given data sets one 
must calculate the detection system sensitivity, (3, given in volts output per unit displacement 
of the particle. This is because the sensitivty can alter between experiments due to variations 
in power or simply geometry at the focus. Finding (3 allows voltage power spectra, those 
directly from data, to be converted to physical spectra, nm 2 Hz _1 versus Hz. Conventional 
methods rely on the relative simplicity of colloidal systems by using, for example, the drag 
force method [49], its extension to an oscillating sample stage [50] or moving a fixed bead 
over a known distance through the laser beam waist [51]. Clearly the former two would be 
difficult to implement in air and the latter is obviously not a good replica of experimental 
conditions [52]. A recent technique has been demonstrated that combines two techniques to 
measure detector calibration from experimentally measured values alone [50]. It is hoped, 
even with the unique problems of airborne trapping, by using AODs or SLMs to oscillate 
the trap position this technique will be developed for future experiments. 

Here we are not concerned with high precision and for simplicity we calculate the detector 
sensitivity, /3, from an uncalibrated voltage power spectrum S v (u) = f3 2 Si nert {uo) using the 
plateau value, P v , reached for uj 3> VL in the function uj i S v {uj) [11]. We find the detector 
sensitivity, (3, to be; 
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III. EXPERIMENTAL 



Droplets are trapped using a custom built inverted tweezers pictured in fig. 1. The 
beam from a 532nm Laser Quantum Finesse 4W c.w. laser is expanded by a Keplerian 
telescope to slightly overfill [53] the back aperture of a Nikon Plan lOOx (NA = 1.25 [54]) 
oil immersion microscope objective. The beam is focused though a type one cover slip 
into an aerosol chamber constructed from a cylindrical plastic enclosure 9mm in height and 
35mm in diameter. This produces an enclosed environment where a high relative humidity 
can exist and also shields the trapping region from external air currents. The top of the 
chamber is made from a type zero cover slip to allow for transmission and then collection of 
the scattered trapping laser by a long working distance (LWD) Mitutoyo lOOx (NA = 0.55) 
objective, whose back aperture is imaged [9] equally onto the four quadrants of a quadrant 
photodiode (QPD) (Hamamatsu Silicon Diode Array S5980) via a 4f imaging system. The 
Mitutoyo objective also acts as the condenser lens for K older illumination (not shown). The 
Nikon objective and an appropriate tube lens images the sample through a laser filter onto 
a Basler A602f firewire camera. 

The liquid aerosol is produced by nebulising a salt solution (20-80g/L) with an Omron 
MicroAir NE-U22 vibrating mesh nebuliser which produces a polydisperse sample of liquid 
droplets with a mass median aerodynamic diameter of 4.9/xm [55]. The aerosol is transferred 
through a hole in the chamber side via a custom made tapered glass nozzle [56]. 

The trapping beam is focused ~ 30yum above the coverslip which are soaked in a 50% 
aqueous dilution of 'Decon 90' for longer than one week. This treatment increases the 
hydrophilicity of the glass and once aerosol has deposited on the cover slip it provides a 
relatively thin, flat, and uniform film of water above which we trap. De-ionised water 
saturated tissue paper is also placed in the chamber to increase the relative humidity, but 
we ensure it does not touch the cover slip as this can induce flows in the water layer. 
Figure 2 shows an enlarged view of the trapping region's geometry and also explains the 
relation between trapping height, L, and objective displacement, X. 

Control over droplet size was required to fully investigate observed phenomena. Firstly, 
this was achieved imprecisely by varying the concentration of the nebulised salt solution [58] 
as a higher concentration decreases the droplets vapour pressure allowing them to equilibrate 
with their surroundings at larger sizes. Secondly, more precise size selectivity can be induced 
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FIG. 1: Apparatus diagram. A Gaussian beam is expanded by lenses LI and L2 and directed to 
slightly overfill the back aperture of the Nikon objective (TOBJ) with mirrors M and DM. The 
long working distance objective (LWD) collects the droplet 's scattered light and its back aperture 
is imaged onto the QPD via a 4f lens system. Power is controlled using a polarising beam cube 
(PBC) and half wave plate (WP). The same Nikon objective with an appropriate tube lens (TL) 
is used to image the sample (S) through a dichroic mirror (DM) and filter (F) onto the firewire 
camera (CMOS). The QPD, COB J, and TOBJ are each mounted on three axis translation stages 
with TOBJ's axial axis controlled either manually or by digital micrometer. BD is a beam dump. 




FIG. 2: Enlarged view of the trapping region in the sample of figure 1. The refractive indices of the 
coverslip and index matched oil, water, and clir £1X6 i^w 

and n a respectively. Ah is the thickness 
of the water layer. Displacing the objective X microns from being focussed on the first interface 
displaces the particle a distance L, given in the figure. 
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with, on average, a positive linear dependence of captured droplet size on laser power [30, 31]. 

Having trapped a droplet the nebuliser is turned off. Once the droplet has reached 
equilibrium with its surrounding environment, and the remaining aerosol settled, the current 
produced by the detection of light on the QPD is sent, via shielded cables, to amplification 
electronics [57] containing a 50 kHz anti-aliasing filter. Data was acquired at a sampling 
frequency of 50 kHz for four seconds with a National Instruments PCI-6014E DAQ card, 
in differential mode. The voltage difference between left and right pairs of quadrants on 
the QPD represents the x position and the difference between the top and bottom pairs 
represents the y position. The voltage versus time data was Fourier transformed using 
Lab VIEW and all remaining data analysis was performed offline at a later time. In order 
to minimise any parameter variation over time the experiments were carried out as quickly 
as possible with raw voltage versus time data not saved to increase speed still further. The 
detailed analysis of the data obtained for a colloidal case is extensively described in Berg- 
S0rensen and Flyvbjerg [37], and much remains the same here. An image of the trapped 
droplet was also taken with each power spectrum for later analysis. 

To further reduce background noise work was always carried out, where possible, solitarily 
in the laboratory with the laser used at > 30% capacity and power control achieved by 
using a pair of half wave plates with polarizing beam cubes. The first split the beam for 
two different experiments and the second controlled power for this experiment alone. The 
power was varied between a minimum of 0.702 ± 0.009mW and a maximum of 510 ± 6mW. 

Unlike tweezing in water, simply increasing the trapping power does not assist in cap- 
turing an aerosol droplet from the nebulised cloud and as such initial laser power must be 
carefully selected. Each droplet trapped was subject to an increase in laser power in uniform 
steps with power spectra measurements taken at each. The minimum attainable damping 
ratio for each droplet was taken from the last power spectra measured before it fell from the 
trap upon increasing the power (i.e. the highest power). This represents an upper limit on 
the ratio for that size. 

To study how the water-air interface to droplet height may affect the dynamics we simply 
keep a constant laser power and vary the height of the sample stage, controlled and measured 
by micrometer. The water layer thickness was measured by observing when a reflection 
of the trapping beam focus is obtained at both the water-air and glass-water interfaces. 
Having been focussed through two index mismatched interfaces (glass to water and water 
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to air) there will be an associated focal shift [59] of which a rigorous description of this is 
complex [60, 61] and not discussed here. A simple paraxial approximation is used to calculate 
the droplets position inside the chamber given a vertical displacement of the sample stage 
around a fixed objective. Modelling of the trap's axial equilibrium position and experiments 
imaging the droplet from the side indicate that the relationship between droplet height and 
objective displacement is linear (data not shown) supporting the paraxial assumption [29, 
63]. 

IV. RESULTS 

Typical power spectra of position fluctuations from optically trapped droplets are shown 
in figure 3, illustrating, for a 3.7 ± 0.2 /im radius droplet, the ease with which the system 
can be transferred between over- and under-damped dynamics by varying laser power. The 
tail falls off with oj~ a as expected for u 3> f2 from equation 2 and a clear resonance peak 
begins to establish itself, indicative of the droplet moving through the critical and into the 
under-damped regime. 

For completeness, a plot of the autocorrelation function [14] of a single droplet in an 
under- and over-damped state is shown in figure 4. It shows the classic exponential decay 
for over-damped motion and sinusoidal oscillation enveloped by exponential decay for under- 
damped oscillators as expected. 

The trend seen in figure 3 remains for all droplets; an increase in power increases lateral 
trap stiffness and moves the system towards or into the under-damped regime. A range of 
damping ratios has been observed from 3.57±0.07 down to 0.260±0.006 over the 4.7±0.5 /xm 
radius range studied. There is also an associated decrease in area under the power spectrum 
curve with increasing laser power, indicating a reduction in the position variance of the 
droplet. 

The inclusion of inertial terms in the Brownian theory means only the particle's mass is 
needed to calculate trap stiffness. Using the radius from video microscopy we obtain lateral 
trap stiffness values ranging from 0.12±0.10 to 98± 17 pN/zm -1 for 1.0±0.3 to 5.7±0.4 /iin 
radius droplets. 

One would expect the natural frequency of trapped droplets to vary as the square root 
of laser power, assuming the trap stiffness is linearly proportional to trapping power. This 
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FIG. 3: Power spectra of a droplet of radius 3.7 ± 0.2 /zm trapped at powers 40.9 ± 0.5 (circles), 
130±2 (squares), and 356±4 (triangles) mW resulting in damping ratios of 1.69±0.04, 0.794±0.01, 
and 0.364 ±0.001 respectively. As the power increases the appearance of a resonance peak is clear 
indicating the move into an under-damped regime, along with a decrease in area and hence position 
variance. The natural frequency increases with laser power because of the associated increase in 
lateral trap stiffness, k. 
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FIG. 4: Plot of experimental autocorrelation data with associated fits for a 5.2 ± 0.2 /im optically 
trapped aerosol in an over- (triangles) and under-damped (circles) state trapped with powers 74 ±1 
and 442 ± 7 mW respectively. Both traces clearly follow the classic exponential decay except in 
the under-damped case there is also the sinusoidal oscillation expected. 

is confirmed in figure 5 for a 1.8 ± 0.2 /im radius droplet. For the range of radii and powers 
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studied here I observe natural frequencies between 27r(328 ± 12) Hz and 2-7r(3433 ± 15) Hz, 
falling close to and well above the corner frequencies measured by tweezers in liquid based 
systems, although obviously not directly comparable. 
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FIG. 5: An example of how natural frequency for a 1.8 ± 0.2 //m radius droplet increases with the 
square root of laser power as expected from Q = \J n/m. The lateral trap stiffness axis is displayed 
for interest and is non linear. The error bars are standard error of the mean for the natural 
frequency rather than the trap stiffness (although they are smaller than the points themselves). 

As mentioned a curiosity of airborne trapping is an upper limit on laser power above 
which particles fall from their trap. To ascertain the cause of this clear instability we plot, 
in figure 6, the upper limit on the minimum damping ratio attainable against droplet radius. 

Figure 6 illustrates there is no problem in transferring from over- to under-damped dy- 
namics with no instability induced, suggesting this is not the reason for the odd phenomena 
seen. There is clear size dependence but one must be careful to note this does not lead to 
the conclusion that small droplets cannot exist in the under-damped regime; this cannot be 
excluded for certain as they fall from the optical traps before they reach it. 

The above results produce a downward shift in damping ratio by increasing lateral trap 
stiffness with larger laser powers. A decrease in friction felt by the droplet could likewise 
shift the ratio by varying the damping and as stated earlier Faxen's correction predicts that 
the proximity of a surface to our microscopic object heavily influences this. Utilising this 
surface to droplet height dependence figure 7 demonstrates that lowering the sample stage, 
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FIG. 6: Upper limit on the minimum attainable damping ratio against droplet radius. It is an upper 
limit as we increase the laser power in finite increments. The dashed horizontal line represents a 
critically damped system. The error bars are standard error of the mean. 

hence increasing the distance, reduces the damping and transfers the system from over- to 
under-damped. Note the resonance peak remains approximately at the same frequency for 
each spectrum as only the damping is changing, not the trap stiffness, contrary to figure 3. 

Extracting damping values from data similar to fig. 7 we can plot the dependence of 
friction upon droplet-surface height to obtain fig. 8 [62]. Here the micrometer raised the 
sample stage in increments of 1/im, decreasing to 0.5/im as the surface was approached. 

In fig. 9 we plot the natural frequency as a function of height from the water layer. There 
is a steady fall off with distance indicating the spherical aberration induced is degrading the 
trap stiffness. Unlike the data of Vermeulen at al. [52] the fit is approximately linear as the 
trap stiffness is independent of the viscous damping. 

V. DISCUSSION 

The reader may notice significantly more scatter in this investigation's results compared 
to experiments using similar techniques in liquid. The reasons will now be discussed briefly 
showing the complex nature of the experiment and the engineering challenges faced to im- 
prove future precision. 

The majority of previous work using the power spectrum method is based on tweez- 
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FIG. 7: Power spectra demonstrating changes in damping as a function of surface-droplet distance 
for a 3.8 ± 0.2 fim radius droplet trapped with 46.3 ± 0.6 mW. The droplet was moved to heights 
of 4± 1 //m (circles), 9± 1 /mi (squares), and 14 ± 1 p (triangles) above the water layer resulting 
in damping ratios of 3.40 ± 0.06, 1.06 =L 0.01, and 0.92 ± 0.01 respectively. The middle and top 
spectra are multiplied by 25 and 200 respectively to displace the data on the y-axis for clarity. 
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FIG. 8: Variation of damping experienced by a droplet, trapped with 6.2 ± O.lmW, as a function 
of surface-droplet height. The droplet is measured from video microscopy to be 2.80 ± 0.2/im in 
radius and found to be 3.28 ± 0.02mn in radius from the fit to experimental damping data. 

ing solid microspheres, with precisely known radii, in the liquid phase allowing very high 
precision studies; indeed, the ability to detect sphere non- uniformity is possible [50]. In 
the studies here a large source of error is measuring the radius using video microscopy 
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FIG. 9: Natural Frequency as a function of distance from water layer for the droplet in fig. 7. 

with the likely errors propagating heavily into some of the systems calculated properties 
(k oc m oc R 3 ). Trapping of solid aerosols, with known radii would remove this problem but 
this is significantly more difficult [56]. In addition, highly precise radius measurements are 
possible via CERS [72] but require sensitive spectrometers. We have shown in a previous 
publication that fitting to the data of fig. 8 with equation 5 allows the droplet radius to be 
measured with good precision [62][ref]. 

It is difficult to determine which individual factor, trap stiffness or damping, contributes 
to the variation in damping ratio for any given experiment. For an individual droplet the 
surrounding conditions can remain relatively constant over the time of a single experiment 
as, with no additional aerosol flow from the nebuliser, the droplet quickly reaches equilibrium 
with its surroundings. To trap another droplet nebulisation must resume where upon the 
chamber conditions can alter. Additional aerosol can settle on the coverslip changing the 
water layer's thickness and hence the optical potential [61, 63, 67] at the trap site together 
with the proximity of the particle to the surface [37, 48, 63]. As mentioned variation in 
trapping power between droplets alters the droplet's height [29] and hence distance from 
the underlying water layer, thus again altering the optical potential and damping. These 
factors contribute to the system's rather complex and difficult analysis. 

The difficult nature (relative to colloidal tweezers) of trapping in air imposes several im- 
portant experimental methods. A long working distance condenser must be used due to the 
aerosol chamber height, but a higher NA lens may have been desirable to improve detector 
sensitivity [68]. Most colloidal experiments use monodisperse suspensions of solid particles 
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thus allowing an arbitrary number of measurement repetitions; often up to 100 power spec- 
tra are averaged. However, with the dynamic system investigated here (the droplets are 
continuously finding an equilibrium with the surrounding environment) the conditions of 
the experiment may not remain constant long enough for repeated measurements to im- 
prove precision, hence the choice of sampling and no averaging over multiple power spectra. 
Also, I am looking at an inherently unstable region with the aim, at times, of losing the 
trapped droplet so, clearly, another particle of the exact same size and composition cannot 
be found. With the current iteration of apparatus there is a clear trade off between speed 
and precision. 

Some studies have used a secondary, independent probe beam to monitor position fluc- 
tuations as this allows greater flexibility and perhaps improved accuracy [69]. We employ 
only a single beam because a very small amount of power is needed to tweeze in air [31] and 
a second beam would significantly alter the potential at the trap site. 

In future studies we suggest that a system including a 'science chamber' be developed 
where many variables can be controlled. A particle could be trapped and transferred to 
such a chamber with relative humidity control, with or without a water layer, and with a 
lower physical profile to enable the use of higher NA condenser optics. Also the mechanical 
stability of our system is not fully optimised so the precision seen could be improved. 

We have neglected any hydrodynamic corrections, specifically effective mass due to en- 
trained fluid [45] and the frequency dependent Stokes' friction [42], and as such it would 
be pertinent to discuss the effects here. Figure 10 plots the ratio of the power spectrum in 
equation 2 to the hydrodynamically correct version in equation 4. Clearly for a given par- 
ticle type and size the correction is significantly smaller when studying aerosols. However, 
the solid line shows the error begins to become significant for aerosols with a radius that 
would be considered relatively large for particles normally used in power spectrum based 
studies in liquid. Should further studies be performed the hydrodynamic correction must be 
investigated to improve accuracy and precision. 

Not discussed in detail here is that Faxen derived his correction for a sphere moving with 
constant velocity which is not the case, so for a complete solution the frequency dependent 
friction should be combined with Faxen's correction [37, 50]. In addition when dealing 
with solid in fluid systems no slipping occurs at the boundary between the two materials 
upon translation. However, the physics involved becomes more complicated when studying 
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FIG. 10: Theoretical plot of S x f S x -h y dro as a function of angular frequency. For a given particle 
radius the hydrodynamic correction is smaller in air (red dot dashed) than in water (blue dashed), 
but, for the large liquid aerosols (black solid) the effect starts to become significant. Trap stiffness, 
k = 2 pN/iin -1 . 

fluid in fluid systems; slip can occur. Due to the possibility of slip at the fluid sphere's 
surface, flow can be induced inside the sphere. This flow causes reduction in the well known 
pre-factor of Stokes' Law [48, 70] and has been taken into account in the data analysis here. 

VI. CONCLUSION 

The parameters governing the dynamics of optically trapped aerosols have successfully 
been described and experimentally investigated. It has been shown that there is no instability 
induced by crossing through the critically damped regime leading to the conclusion that the 
phenomena seen are due to the optical potential of the trapping laser. The investigation has 
provided results extending the boundaries of precise studies of Brownian motion in optical 
tweezers into a new damping regime. 

It is hoped these results will provide researchers with a new understanding of an old tool 
for studies in both fundamental and applied science, providing a rich playground of study 
in the under-damped regime. 

From the Langevin equation it is seen that there are only four processes providing forces 
that give rise to droplet position fluctuations. There is Brownian white noise, friction, 
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inertia, and the optical force. Transferring from over- to under-damped dynamics does not 
inhibit the stability of the system yet there is clear size dependence with which the droplets 
fall from the trap. The logical conclusion is that the optical force must determine whether 
the droplet remains trapped or not and this will be discussed in a future publication. 
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